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The nature of the underdoped pseudogap regime of the high-temperature superconductors has been 
a matter of long-term debate [1— 3j. On quite general grounds, one expects that due to their low 
superfluid densities and short correlation lengths, superconducting fluctuations will be significant for 
transport and thermodynamic properties in this part of the phase diagram|4, 5j. Although there 
is ample experimental evidence for such correlations, there has been disagreement about how high 
in temperature they may persist, their role in the phenomenology of the pseudogap, and their sig- 
nificance for understanding high-temperature superconductivity [61- 410] . In this work we use THz 
time-domain spectroscopy (TTDS) to probe the temporal fluctuations of superconductivity above the 
critical temperature (T c ) in La2 x Sr x Cu04 thin films over a doping range that spans almost the entire 
superconducting dome (x — 0.09 to 0.25). Signatures of the fluctuations persist in the conductivity in 
a comparatively narrow temperature range, at most 16 K above T c . Our measurements show that su- 
perconducting correlations do not make an appreciable contribution to the charge transport anomalies 
of the pseudogap in LSCO at temperatures well above T c . 

In general, continuous phase transitions are typified by fluctuations with correlation length and time scales that 
diverge near T c . Dynamical measurements like TTDS are a sensitive probe of the onset of superconductivity [TT] and 
measure its temporal correlations on the time scales of interest. In the presence of superconducting vortices such high- 
frequency measurements are not affected by effects like vortex pinning, creep, and edge barriers that often complicate 
interpretation of low frequency and DC results. In this study, we investigate the fluctuation superconductivity in thin 
films of LSCO grown by molecular beam epitaxy (MBE). This synthesis technique provides exquisite control of the 
thickness and chemical composition of the films; the intrinsic chemical tunability of LSCO allows us to investigate 
essentially the entire phase diagram. For details on the films, see the 'Supplementary Information' (SI). 

In Figs.JI^, and b, we show the real (<j\) and imaginary (02) parts of the THz conductivity measured at a number of 
different temperatures for optimally doped LSCO (x — 0.16) with T c =41 K. We obtain similar data at other doping 
levels. The spectra are easily understood in the limiting cases of high and low temperatures. Well above the onset of 
superconductivity, the real part of the conductivity is almost frequency independent and the imaginary part is small, 
consistent with the expectation for the behavior of a metal at frequencies well below the normal state scattering rate. 
At the lowest temperature the conductivity is consistent with that expected for a long-range ordered superconductor; 
o~ 1 is small as most of the low frequency spectral weight has condensed into the uj — delta function, and the frequency 
dependence of cr 2 is very close to l/u. Our principal interest, however, is in the interesting transition region around 
T c , where fluctuations of superconductivity are apparent. Here, both components of the conductivity are enhanced, 
with <7i first rising and then falling as spectral weight moves to frequencies below the measurement range in the 
superconducting state. 

The fluctuation regime is perhaps most evident in Figs. [TJj and d, where we plot the complex conductivity of the 
same sample as a function of temperature. Above the transition, the real conductivity shows a slow decrease as 
the temperature is raised, consistent with the increasing normal state DC resistivity. Near T c a 'loss peak' occurs 
due to the onset of strong superconducting fluctuations; it is exhibited at progressively lower temperatures as the 
probing frequency is decreased (inset to Fig. Ilk). As discussed below, this is a direct consequence of the slowing 
down of superconducting fluctuations as T is lowered. About 10 K above T c , the imaginary conductivity shows a 
sharp upturn at low frequency (inset to Fig. [TJi) indicating the onset of strong superconducting correlations 1 1 2 j . 
Whether such correlations persist well above the 10-15 K range above T c is a more subtle issue, which we discuss 
below. For all dopings, the region of obvious enhancement is far below the temperature of the diamagnetism and 
Nernst onset measured in LSCO crystals [3 EH], but is consistent with the onset found in other AC conductivity 

studies on mon. 

An important quantity for understanding superconducting fluctuations is the phase stiffness, which is the energy 
scale for introducing twists in the phase 4> of the complex superconducting order parameter Ae l ^ r \ As in any 
continuous elastic medium one can write the energy of a phase deformation in the form of E = ^70(V0) 2 where 7^ is 
a stiffness constant. Since a phase gradient is associated with a superfluid velocity, this 'elastic' energy is equivalent 
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FIG. 1: THz conductivity of an x — 0.16 La2- a: Sr a: Cu04 thin film, (a), (b) Real and imaginary conductivities as a function of 
frequency at different temperatures for an optimally doped x = 0.16 sample (T c =41 K). (c), (d) Real and imaginary conductivities 
as a function of temperature at different frequencies for the x — 0.16. In panels (a) and (b) the green curve denotes T c . In panels 
(c) and (d) the vertical lines represent T c . Insets to (c) and (d) show expanded views of the fluctuation region. Black dots mark the 
approximate peak value in the inset to (c). 



to the center of mass kinetic energy. For N particles of mass m, the stiffness is 7^ = One can measure this 

quantity directly through the imaginary part of the fluctuation conductivity 03/ , as k^T^ = — g~ ■ Here is the 

2D stiffness of a single CuC>2 plane given in units of degrees Kelvin, Gq — e 2 /H is the quantum of conductance, and 
t is the inter-Cu02 plane spacing. The phase stiffness is usually regarded as an equilibrium quantity; the effect of 
measuring it dynamically at finite frequency uj is to introduce a length scale over which the system is probed. In 
models with vortex proliferation this is typically proportional to the vortex diffusion length y ' D juj where D is a 
vortex diffusion constant. 

In Fig. [2j we plot the quantity cjer 2 as a function of frequency and temperature, (a) and (b) respectively, for 
an underdoped sample (x = 0.09) with T c =22 K. This quantity is proportional to the phase stiffness when the 
superconducting signal dominates over the normal state background. Due to uncertainties in its possible form, we 
have abstained from subtracting a background from the plotted signal. We find, however, that different choices for 
backgrounds have little affect on our conclusions. Deep in the superconducting state, there is essentially no frequency 
dependence to ojct^, which is consistent with the fact that the system's phase is 'stiff' on all length scales. At higher 
temperatures the curves in Fig[2]D spread as a frequency dependence is acquired. Interestingly, the temperature where 
the spreading first becomes significant is very close to the temperature where the Kosterlitz-Thouless-Berenzinskii 
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FIG. 2: uo2 is proportional to the THz 'phase stiffness' (a) UJ02 vs. frequency at different temperatures for the x = 0.09 
underdoped sample (T c =22 K). The green curve denotes T c . The pink curve is the effective KTB temperature, (b) UJ02 as a function 
of temperature for different frequencies for the x = 0.09 underdoped sample. The dotted black line is the prediction of the KTB 
transition for a single isolated CuC>2 plane. The temperature of its crossing with the experimental data define an effective KTB 
temperature. The vertical line represents T c . (c) uj(J2 at 800 GHz at a number of different dopings. T c is indicated by the asterix 
and T is indicated by diamonds, (d) T (red diamond) and error bars (red lines), as defined for the x = 0.16 sample by the onset 
of curvature of cjct2. 



(KTB) theory (black dashed line) would predict a discontinuous jump in this quantity for an isolated CuC>2 plane 
|18j . LSCO is of course only a quasi-2D system, but the fact that we observe a spreading at the KTB prediction 
shows that at some temperature T^tb < -^c there begin to be significant fluctuations of a 2D character even below 
the transition. T c itself does not occur until higher temperatures, as it is controlled by 3D couplings. All of our 
underdoped samples, as well as previous measurements of BSCCO films [11] show this behavior. As evidenced by the 
frequency dependence of ujo-i above T^ B , the fluctuations first degrade the stiffness on long length and time scales, 
i.e. low frequencies. To compare across the phase diagram, in Fig. [2j: we plot w<72 for many dopings at a frequency 
0.8 THz. Qualitatively, we regard the onset temperature T a as the temperature where the quantity lou^ presents a 
substantial deviation from the trend of the higher-temperature normal state. As a quantitative measure, we find that 
all samples exhibit a very sharp deviation out of the smooth high-temperature signal in this range in plots of the 
second derivative (Fig. [2ji) vs. temperature. We take this deviation as the quantitative measure of T D (with error 
bars). For all dopings, T Q is no more than 16 K above T c . It is also interesting to note that the difference between T Q 
and T c is relatively constant over the experimental range. 
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FIG. 3: Characteristic fluctuation frequency obtained through a scaling analysis, (a) Phase of the conductivity ip — tan^o^/ci 

in the frequency range from 0.5 THz - 1.5 THz as a function of the reduced frequency for 45 different temperatures in the range 
of 30 K - 55 K for the x = 0.16 sample, (b) The superconducting fluctuation rate fi(T) determined as detailed in the text for a 
variety of different dopings. Colored vertical lines are the T c 's for the respective dopings. It is observed that the fluctuation rate for 
all dopings approaches a limiting linear temperature dependence. 



As mentioned above, the data in Figs. [TJ and [2] are consistent with superconducting correlations that are typified 
by a slowing down of the characteristic fluctuation rates as temperature is decreased. In general, the diverging length 
and time scales near a continuous phase transition lend themselves to a scaling analysis in which response functions 
can be written in terms of these diverging scales. Close to a superconducting transition one expects that the relation 

holds for the portion of the conductivity due to superconducting fluctuations 07. Here T® is a temperature dependent 
prefactor and f2 is the characteristic fluctuation rate. Temperature dependencies enter only through the quantities £1 
and T?. This scaling function is similar to the one proposed by Fisher, Fisher, and Huse (FFH) [TU] and is identical 
to the one used in previous THz measurements on underdoped BSCCO [11]. Note that Eq. [T] is a very general 
form that does not assume any particular dimensionality of the system or character (vortex, Gaussian, etc.) of the 
fluctuations or functional dependencies on temperature of fi and T°. Also note that while the fluctuation conductivity 
at = \<7f\e %v is a complex quantity all prefactors to the S scaling function are real and therefore the phase of S is 
equal to the phase of o~f. S is expected to exhibit single parameter scaling and thus a collapse of ip measured at 
different temperatures as a function of the reduced frequency uj/fl yields the temperature dependent fl. In Fig. [3^,, 
we show the collapsed phase tp = tan -1 ^/ '<J\ of the x = 0.16 sample as a function of the reduced frequency uj/il for 
45 different temperatures in the 30 K to 55 K range. As expected, the phase is an increasing function of w/fi, with 
the metallic limit ip = reached at ui/fl — > and <p becoming large (but bounded by 7r/2) as ui/fl — > 00. 

We were able to perform the scaling analysis and obtain similarly good data collapse for all eight samples in the 
range x = 0.09 to 0.25. In Fig. [3Jd we plot their extracted fluctuation rates fi as a function of temperature. In 
almost all samples the fluctuation rate increases very quickly (almost exponentially) within 10 K of T c . Note that the 
characteristic fluctuation frequency fl is actually smooth through T c , in a manner quite different than expected from 
the FFH scaling. And it appears to extrapolate to zero near the effective Tfftf B where 2D-like fluctuations become 
strong. At a temperature we denote as Tq, the extracted fluctuation rate crosses over to a regime where it grows at 
a much smaller rate that is proportional to temperature as aksT/h. As the resistivity itself is linear in T (See SI), 
Tq is an alternative measure of the temperature scale where we cannot distinguish superconducting fluctuations from 
normal state transport. Interestingly, we continue to get good scaling and data collapse as we continue the analysis 
for another 5-10 K above Tq. This behavior may be related to the cj/T scaling that has been seen in the normal state 

In Fig.[4]we construct a phase diagram of the fluctuation regime summarizing our results. We compare T c of films and 
bulk crystals, the diamagnetism onset from Ref. [TU] , the region of THz conductivity onset T a (as defined by the sudden 
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FIG. 4: Phase Diagram of the onset of superconducting correlations in LSCO. Temperature vs. doping phase diagram of 
La2-iSr^Cu04 comparing T c of thin films and bulk crystals with the THz conductivity onset T and the diamagnetism onset 
temperature from Ref. |10| . Here, T a is expressed as a shaded region to convey the uncertainty in its determination. The temperature 
Tq at which the characteristic fluctuation rate Q becomes proportional to temperature is plotted in green. 



change in curvature of oj<72 from Fig. j2p), and Tq where Q begins to grow linearly with temperature. Unsurprisingly, 
to within experimental uncertainty Tq and T a track each other . For all dopings the THz fluctuation conductivity 
shows an onset between 5 K and 16 K above T c . This is in strong contrast to measurements like diamagnetism, 
where the signal persists in similar samples to almost 100 K above T c . Although fluctuation diamagnetism is a 
thermodynamic quantity sensitive to spatial correlations and THz conductivity a dynamical quantity sensitive to 
temporal correlations, the differences are still surprising as within most models based on diffusive dynamics, one 
expects a close correspondence between them [18 . An approximately 20 K difference in the onset temperatures exists 
between diamagnetism and previous THz measurements on BSCCO dULU], but the present 100 K discrepancy is 
much larger (particularly when normalized to the respective T c 's) and cannot be easily explained away. 

We may draw a few possible conclusions. One obvious possibility is that there is little correspondence of diamag- 
netism (and Nernst) with conductivity because the former class of measurements are sensitive to something other than 
only superconductivity well above T c . In this regard, it has been shown recently that the onset of density wave order 
can give a strong Nernst response [9], and that bond currents states can in principle exhibit enhanced diamagnetism 
[22) . However, if the diamagnetism signal above T is solely due to superconducting correlations then it is a well-posed 
theoretical challenge to explain the lack of straightforward correspondence with conductivity. 

We have found that at some temperature Tq > T c the fluctuation rate Q is either overwhelmed by the linear-in-T 
normal state scattering rate or becomes linear in T itself. Calculations that model the normal state as a vortex 
liquid [23H25) with a characteristic dissipation rate proportional to T favor the latter scenario. Unfortunately our 
measurements cannot distinguish these possibilities. It is interesting to note that the regime in which we observe a 
large fluctuation conductivity is essentially the same as the regime of "fragile London rigidity" in the diamagnetism 
of Li et al. [T0|. I21j. It may be that there are two distinct types of superconducting fluctuations; one of which gives 
a divergent or near divergent contribution to the susceptibility in a relatively narrow range of temperatures above 
T c and the other of which gives a more extended, but less spectacular contribution to the susceptibility. If the first 
type makes a more significant contribution to the conductivity than the second, this may resolve the apparent conflict 
between the two probes. Similarly, the differences may arise in how the quantities in question depend on correlations 
in space (probed by diamagnetism) vs. correlations in time (probed by conductivity). As noted above, within classical 
diffusive dynamics, one generally expects a correspondence between these quantities. Our data show that if the regime 
of enhanced diamagnetism is due to superconductivity then a very unconventional relationship must exist between 
length and time correlations in the fluctuation superconductivity. Among other possibilities, this could arise from 
phase separation [55], unusually fast vortices [S7J HB], or the presence of explicitly quantum diffusion. Whatever 
the reason, our data show that superconducting fluctuations do not make an appreciable contribution to the charge 
transport anomalies in the pseudogap regime at temperatures well above T c . 
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I. METHODS 

The complex conductivity was determined by time-domain THz spectroscopy. A femtosecond laser pulse is split 
along two paths and excites a pair of photoconductive 'Auston'-switch antennae grown on radiation damaged silicon 
on sapphire wafers. A broadband THz range pulse is emitted by one antenna, transmitted through the LSCO film, 
and measured at the other antenna. By varying the length-difference of the two paths, we map out the entire electric 
field of the transmitted pulse as a function of time. Comparing the Fourier transform of the transmission through 
LSCO to that of a reference resolves the full complex transmission. We then invert the transmission to obtain the 
complex conductivity via the standard formula for thin films on a substrate: T(w) = i +n +Zoa{ui)d e ' 1 ^''' wnere is 
the phase accumulated from the small difference in thickness between the sample and reference substrates and n is 
the substrate index of refraction. By measuring both the magnitude and phase of the transmission, this inversion to 
conductivity is done directly and does not require Kramers-Kronig transformation. 

Our scaling analysis includes some uncertainty in setting the overall scale of ft as Eq. [I] only specifies ft up to a 
numerical factor. We set the overall scale of f2 so that the loss peak in Fig. [Tp is exhibited at a temperature when 
the probing frequency u> equals the fluctuation frequency at that temperature. This is an imprecise procedure 
and applying it strictly gives a small distribution (40%) in the values of a for different dopings. We choose the 
normalizations for the data plotted in Fig. [3] such that a = 0.4 which is the mean value for all the data. We take the 
relatively small spread in a's as evidence of the veracity of our procedure but emphasize that the ct's could be revised 
with a different criteria for the scale factors. As a crosscheck we also fit the conductivity using a Kramers-Kronig 
consistent fitting routine and get good agreement above T c between the half width of the Lorentzian peak used 
to model the fluctuation contribution and the rates given in Fig. [3]d. 

The LSCO films were deposited on 1-mm-thick single-crystal LaSrA104 substrates, epitaxially polished perpendic- 
ular to the (OOf) direction, by atomic-layer-by-layer molecular-beam-epitaxy (ALL-MBE) [? ]. The samples were 
characterized by reflection high-energy electron diffraction, atomic force microscopy, X-ray diffraction, and resistiv- 
ity and magnetization measurements, all of which indicate excellent film quality. For accurate determination of the 
conductivity, it is critical to know the film thickness accurately. This was measured digitally by counting atomic 
layers and RHEED oscillations, as well as from so-called Kiessig fringes in small-angle X-ray reflectance and from 
finite thickness oscillations observed in XRD patterns. The x=0.f4, 0.16, and 0.25 films are 80 monolayers thick; the 
x=0.095 and 0.19 Aims are 114 monolayers; the x=0.09, 0.10 and 0.12 films are 150 monolayers (one monolayer is 
» 6.6A ). 
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IV. SUPPLEMENTARY INFORMATION: FILM CHARACTERIZATION BY XRD, AFM, MUTUAL 

INDUCTANCE AND RESISTIVITY MEASUREMENTS 

We have characterized all samples by X-ray diffraction and low-angle X-ray reflectance using an X-pert PRO 
PANalytical four-circle diffractometer with a Cu tube. A typical X-ray diffractogram for the x = 0.1 sample is 
shown in Figure [5] Out-of-plane lattice constants Co are determined using multiple Bragg diffractions and applying 
the Nelson-Riley algorithm (1). For the various samples, the nominal thickness is known by virtue of our digital 
(atomic-layer-by-layer) molecular beam epitaxy technique. This nominal thickness was in agreement to within ±5% 
with the thickness calculated from the finite-thickness oscillations (so-called Kiesig fringes) seen in the diffractograms 
as satellite peaks on both sides of the strong Bragg peaks (see the inset to Figure [5| , as well as in low-angle X-ray 
reflectance. The error in thickness determination (typically of the order of 5-10 A) depends on the accuracy that one 
can measure the period of oscillations. The large number of well-resolved satellite peaks in Fig. [5] indicates that the 
surface roughness of the film is small. 
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FIG. 5: X-ray diffractogram of the LSCO x = 0.1 film. An accurate value of the lattice constant Co can be retrieved by means of 
the Nelson-Riley algorithm using a number of observed diffraction peaks. Inset: Magnified view of the (004) Bragg peak, showing 
pronounced finite-thickness oscillations. 

This small surface roughness was confirmed with Atomic Force Microscopy (AFM). For the 10 x 10 /xm 2 field of 
view shown in Fig. [6] for the x = 0.1 sample, the typical rms surface roughness (2) varied between 0.2 and 0.9 nm; in 
all cases it was less than one unit cell height (en ~ 1.3 nm). 

Every sample was also characterized using a two-coil transmission mutual inductance (MI) technique to check 
film homogeneity. The diamagnetic MI signal onsets when the resistance of the sample reaches zero. The in-phase 
(dissipative) MI shows a peak below T c , the width of which is sensitive to both in-plane and out-of-plane inhomogeneity 
up to the scale of the coil diameter (w 1 mm) . As seen in Fig. [7^, for a optimally doped sample the measured FWHM 
values for the peaks were in the range between 0.2-0.5 K, which demonstrates the high homogeneity of the films and 
their excellent quality. Any variation in stoichiometry, either in-plane or out-of-plane, would result in broadening of 
the dissipative peak, or even in the appearance of multiple peaks. 

Measurements of resistance as a function of temperature provide another figure of merit attesting to the high 
quality of the films, the Residual Resistivity Ratio (RRR). As seen in Fig. [7J;, for an optimally doped LSCO film the 
resistivity p(T) is typically linear in T down to temperatures near T c . This linear dependence can be extrapolated 
down to T = 0, to define the so-called Residual Resistivity p - The ratio p^/ p(T=300 K) can be used as a measure 
of film quality. In general a RRR of 5-10 % is considered very good. In this particular sample the RRR <C 1%. 
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FIG. 6: Typical AFM image of the z = 0.1 LSCO films. The surface of this sample show terraces that are one mono-layer (1/2 
unit cell) high. The rms roughness extracted from this 10 x 10 fim image is 2.74 A 
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FIG. 7: Transport in LSCO thin films, (a) The in-phase (dissipative) and (b) the out-of-phase (inductive) components of mutual 
inductance measured by two-coil technique, (c) Resistivity vs. temperature measured in an optimally doped LSCO film. The red 
straight line is the linear fit from the room temperature down to few degrees above T c . 
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